Astronomy & Astrophysics manuscript no. deflection' v3 


©ESO 2010 


August 18,2010 





o 

(N 
< 
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Multiple-deflection, reduced shear and lensing bias corrections 
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ABSTRACT 

It is usually assumed that the ellipticity power spectrum measured in weak lensing observations can be expressed as an integral over 
the underlying matter power spectrum. This is true at order (9(0 2 ) in the gravitational potential. We extend the standard calculation, 
constructing all corrections to order (9(<I> 4 ). There are four types of corrections: corrections to the lensing shear due to multiple- 
deflections; corrections due to the fact that shape distortions probe the reduced shear y/(l - k) rather than the shear itself; corrections 
associated with the non-linear conversion of reduced shear to mean ellipticity; and corrections due to the fact that observational galaxy 
selection and shear measurement is based on galaxy brightnesses and sizes which have been (de)magnified by lensing. We show how 
the previously considered corrections to the shear power spectrum correspond to terms in our analysis, and highlight new terms that 
were not previously identified. All correction terms are given explicitly as integrals over the matter power spectrum, bispectrum, and 
trispectrum, and are numerically evaluated for the case of sources at z = 1. We find agreement with previous works for the (9(0 3 ) 
terms. We find that for ambitious future surveys, the (9(0 4 ) terms affect the power spectrum at the ~ 1 - 5cr level; they will thus need 
to be accounted for, but are unlikely to represent a serious difficulty for weak lensing as a cosmological probe. 

Key words, cosmology: theory - gravitational lensing - large scale structure of the Universe 

1. Introduction 
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jH ' Cosmic shear, the distortion of light from distant galaxies by the tidal gravitational field of the intervening large scale structure, is 
c/2 , an excellent tool to probe the matter distribution in the universe. The statistics of the image distortions are related to the statistical 
, ^ , 1 properties of the large scale matter distribution and can thereby be used to constrain cosmology. Current results alread y demonstrate 
the power of cosmic shear observations at constrain i ng the clusteri ng amplitude <x 8 and the matter density Q m (e.g.. lMassev et all 
{SJ . l2007bllSc hrabback et al.l2007tlBeniamin et al.l2007trFu et al.l l2008). Furthermore, cosmic shear provides an ideal tool to study dark 
' energy through measuring the evolution of non-linear structure with large future surveys (DESQ LSSTBJDEIvQ, EuclicQ). These 
^sO \ upcoming large weak lensing experiments will limit the statistical uncertainties to the percent level. 

In order to extract cosmological information from these cosmic shear experiments, the increased data quality needs to be 
! ' accompanied by a thorough treatment o f systematic err ors. On the observational side, this requires accurate information on the 
^] , redshift distribution of source galaxies dMa et al.ll2006l) and precise measurements of galaxy shapes which correc t for observa- 
■ tional systematics such a s pixelization, noise, blurring by seeing and a spatially variable point spread fun ction (se e Massev et al.l 
y— I ' | 2007at iBridle et alj|2009l) . On the theoretic al side, astrophysical co ntaminants, like source lens clustering (Bernard eau et al.lll997l: 
0^ . ISchneider et al.l2002l) , i ntrinsic alignmen t ( King & Schneiderf2003l) and the correlation between the gravitational shear and intrinsi c 



■ ellipticities of galaxies dHirata & SeliakH2004t lKingll2005l;lJoachimi & Schneiderfl2008l:IZhangl l2008: Joachim ilTSchneideif2 009). 



\ need to be understood and removed. The prediction of lensing observables also requires precise models of the non-linear matter 
• *h , power spectrum and models for the relation between lensing distortion and large scale matter distribution which go beyond linear 
/\ ' theory. While N-body simulations m ay predict the non-linear dark matter power spectrum with percent level accuracy in the near 
^ | future dHeitmann et alj 2008L 20091). the eff ect of baryons, which is a significant contamination to the weak lensing signal above 



I ~ 2000 (Jing et al. 200d lRudd et aTll2 0081. is more difficult to account for and is the subject of ongoing work. 



In this paper, we consider corrections to the relation between the observed lensing power spectra and the non-linear matter 
density field. In the regime of weak lensing, the observed galaxy ellipticities (<?/) are an estimator of the reduced shear gj = Ji/(1-k), 

(e,) = C-^- , (1) 

1 - K 

where C is a constant which depends on the type of ellipticity estimator (e.g. ISchneider & Seitzlll99"5t ISeitz & Schneiderfll997l) 
and the properties of the galaxy population under consideration, ji is a component of the shear, k is th e convergence, and the 
subscript / refers to the two components of the ellipticity/shear (see e.g. iBartelman n & Schneidedl2001l for more details). The 

1 www.darkenergysurvey.org 

2 www.lsst.org 

3 http://jdem.gsfc.nasa.gov 

4 http://sci.esa.int/science-e/www/object/index.cfm?fobjectid=42266 
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two-point statistics of the measured ellipticities are simply related to the reduced shear power spectrum. ICoorav & Hul ((2002) have 
calculated the shear power spectrum to fourth order in the gravit ational potential. For the reduced she ar power spectrum there 
exists an approximation to third order in the gravitational potential (Dod elson et al . 2006). Shapiro (2009) has demonstrated that on 
angular scales relevant for dark energy parameter estimates the difference between shear and reduced shear power spectra is at the 
per cent level and ignorin g these corrections will noticeably bias dark energy parameters inferred from future weak lensing surveys. 

ISchmidt et alj d2009al) introduced another type of corrections, termed lensing bias, which has a comparable effect on the shear 
power spectrum as the reduced shear correction: Observationally, shear is only estimated from those galaxies which are bright 
enough and large enough to be identified and to measure their shape. This introduces c uts based on observed brightness and observed 
size, both of which are (de)magnified by lensing (e.g. Broad hurst et al.lfl995b lJainll2002h . and will thus bias the sampling of the 
cosmic shear field. 

In the following we complete the calculation of the reduced shear power spectrum to fourth order in the gravitational potential 
to include multiple deflections and to account for the effects of lensing bias and the non-linear conversion between ellipticity and 
reduced shear. We consider all lensing-related effects through (9(<t> 4 ), but do not include effects associated with the sources (source 
clustering and intrinsic alignment corrections). 

This paper is organized as follows: We describe our technique for calculating higher order lensing distortions and power spectra 
in Sect. [2] Derivations of the different types of corrections to the shear and reduced shear power spectra are given in Sect. l3.1l through 
Sect. 13.41 We quantify the impact of these corrections on future surveys in Sect. [4] and discuss our results in Sect. [5] 



2. Calculational method 

In this section we derive the higher order lensing distortions following iHira ta & Seljak| (l2003h . and introduce our technique and 
notation for calculating power spectrum corrections. 

Throughout this calculation we assume a flat universe and work in the flat sky approximation. We use a unit system based on 
setting the speed of light c = 1, which makes potentials dimensionless. We use the Einstein summation convention and sum over 
all Roman indices appearing twice in a term. Lower case, italic type Roman indices a,b,c, ... = 1,2 are used to for Cartesian 
components of two dimensional vectors and tensors; capital case, italic type Roman indices I,J,K, ... = 1,2 are used for the 
components of polars which are defined with reference to a Cartesian coordinate system but have different transformation properties. 
Greek indices are used forredshift slices. 



2. 1 . Lensing distortion tensor 

We work in the flat sky approximation and choose the sky to lie in the xy-plane. Photons travel roughly along the — z direction and 
are deflected by the Newtonian potential <£ generated by the nonrelativisti c matter inhomogeneities. As long as their deflection from 
the -z direction is small, they observe a metric (e.g. lHirata~& Seliak 2003) 

ds 2 = a 2 (r) [-(1 + 2<&)dr 2 + (1 - 2<D) [d X 2 + x 2 (dn 2 x + dn 2 ))] , (2) 

where a is the scale factor, x is the comoving radial distance, and n is the angular coordinate of the photon path on the sky. We 
calculate the deflection angle of a light ray from its null geodesic equation 



_d_/dn \ 

d X U/j 



^ d®(x(n,x);z(x)) m 
- 2 dn * ' (3) 



where <5(x; z) is the Newtonian potential at position x and redshift z, with initial conditions n(x = 0) = no and d x n(x — 0) = 0. 

To first order in cD, the integration is performed along the unperturbed photon trajectory, this is the so-called Born approximation. 
Taylor expanding Eq. (f3]l to third order in <t> we obtain a perturbative solution for the deflection angle d = n - no 

m(z s ) = n 0; + flffc) + df\ Zs ) + df\ Zs ) 

=n 0i + d m (z s ) - 2 £dxW(x,Xs)x 2 ®,ia(x)d ( V(x) - 2 £d x W(x,Xs)x 2 (^aaCtf^Ctf^fr) + *»0r)<g°(*)) (4) 

= «o, -2 dxW(x,Xs)X®i(X)+4\ d x W(x,x s )x 2 \ dx'W(x',x)x'®j a (x)®,a(x') 
Jo Jo Jo 

4 d X W(x,x,)x 3 dx'W(x',x)x' d X " W( X ",X)x"®iab(x)®a(x')W) 
Jo Jo Jo 

+ 8 dxW(x,X,)x 2 dx'W(x',x)x' 2 d X " W(x" ,x')x''®ia(x)®ab(x')W) . (5) 
Jo Jo Jo / 

where W(x',x) — ("t _ ~) ®(x ~ X') w i m ®00 the Heaviside step function. Here X s = X(Zs) 1S the comoving distance of a source 
at redshift z s , commata represent comoving spatial transverse derivatives. These spatial derivatives are evaluated at the unperturbed 
position <3>(x) = O(no^,^; z(x)) unless otherwise indicated. The first and second order deflection angles are identical to those found 
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bv lHirata & S eliak (2003)0. The third order deflection angles are caused by the two types of second order transverse displacement 
in the Taylor expansion of <t>(x; z) shown in Eq. (@). We discuss the difference between these terms after Eq. ([8]). 
The distortion of a light ray is then described by the Jacobian matrix 

A(no,z s ) = — — = ' i ,._■_„,. . (6) 



<9n \ ~72 + ^ 1 - k + 7i 

where 7/ are the cartesian components of the shear, and to induces an (unobservable) rotation of the image. Using (0, the distortion 
tensor (fry = 6/j - Ay is given by 



where 



(fry(n , z s ) = ^'(no, z s ) + (frjJW z s ) + (fr^fW Zs ) + (frg B) (n , z.) + ^fW z s ), (7) 
^JWz.) = 2 J o (frWCv,*.)* 2 *^) , 

^, 3A) (n ,z s ) = +4 d^WCr,^)/ <fcW,*)*' W W(x",x)x" ®ijab(x)® a (x')W) , 
Jo Jo Jo 

=>x'<x, x"<x 
r*Xs r*x r*x 

^fWs) = + 8 d X W( X ,x s )x 3 dx'W( X ',x)x' 2 d X " W(x", X )x"^iab(x)^aj(x')W) , and 
Jo Jo Jo 

=^x'<x,x"<x 

r*x& ex ex' 

^■ C) (n ,z s ) = + 8 d X W(x,x*)x 2 dx'W( X ',x)x' 2 d X " Wiy'^x"^— [®ia(x)®ab(y)W)] , (8) 
; Jo Jo Jo on j 

=*x"<x'<x 

where we have used the symmetry of the integrals over^-' and^" in the de rivation of (fr; ; . This calculation automatically includes 
the "Born correction" and "lens-lens coupling" corrections considered by ICoorav & Hul (120021) . Compared to their approach, we 
find additional terms i/^ 3c) which give the third order corrections caused by three lenses placed at different locations along the line 
of sight jx" < y' < y), namely the derivatives of the last term in Eq. ([8). These include the two terms previously considered by 
Shap iro & Cooravl(l2006l) . however, we will show in Sect. B.ll that within the Limber approximation, the 3C term does not contribute 
to the shear power spectrum at (9(<D> 4 ). 

The convergence, shear, and rotation are expressible in terms of (fry by the usual rules k = l(iffu + ^22), 7i = 5(^11 _ tyzi), 

72 = 2"(tAi2 + ^21), and w = ~((fri 2 - 0"2i)- 

Note that while our derivation of the deflection angle is based on the small angle approximation d <sc 1, in the flat sky approxi- 
mation the elements of the distortion matrix need not be as small. 



2.2. Fourier space: first order 

Since we work in terms of power spectra, we need to transform these equations to Fourier space. In the flat-sky approximation, 

(fry(n 0> z s ) = f ^(A l7 (l,z s )e iln ». (9) 
J (2tt) / 

The angular cross power spectra of two fields F and F' is then defined by (F(l)f '(1')) = (27r) 2 Crr'(0f>D(l + 1') with 6d the Dirac delta 
function, which ha s units [<5p(x)l = f xl ~" wh ere n is the dimension of x. Potentials are functions of a three dimensional position 
variable. Following iDodelson & ZhangTd2005l ). we use <p to denote the Fourier transform of the potential in the angular (transverse) 
variables only 

4> d;x) = ^2 / fe 3;zCkr))e^ . (10) 

Then the spatial derivatives of the potential can be expressed in terms of the angular Fourier transform as 

\ M f d 2 l 

® A k...i M (nox;x) = -m 7T - 2 hA---k M ~^x)^ no - (ID 



5 O u r notation differs from lHirata~& Seliak (2003) in using spatial instead of angular derivatives to simplify comparison with Cooray & Hu 
d2002h : IDodelson et alJ l l200lih : IShapird i2tM) 
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Applying this to the first term from Eq. ([8]) and using the relation between convergence, shear and iftjj, we arrive at the well-known 
first order results for convergence and shear 

* (1) fl, z.) = 2 Wd. *) + ^d, z.)) = -/ 2 J o d^WCr,A- s )0(l;^) and yffl, Zl ) = 77(1)* (1) (1, z.) . (12) 

Here Ti(l) = cos(20/) and ^(l) = sin(20;), where tpi is the azimuthal angle of 1. 

We generally decompose the shear components into tangential (or E-mode) shear, je and cross (or B-mode) shear, ys, 

Ted, Zs) = SuTjQtfjO, z s ); y B (l, z s ) = 6/7^(1)^(1, Zs ) , (13) 



with e/y the two dimensional Levi-Civita tensor. To first order, y^ (l, z s ) _ £ ( ^(1, z s ) andy^ (l, z s ) = d/ '(1, z s ) - 0. Their power 
spectra can be obtained under the Limber approximation (Kaiser 1992; Dodelson & Zhang 2005, Eq. (15)), 



(4>(l;xW;x')) = (2nfS D (\ + Y) Sd(X ^ P 9 Q/ X ;z(x)) , (14) 

X 

where P® (l/x', z(x)) is the three dimensional power spectrum of the potential at redshift z(x)- The lensing tomography cross spectra 
between two source redshift slices at z a and zp (with z a < Zp) then read 

Cg l \l;z a ,zp) = C^(l;z a ,zp) = I 4 jT *X W(x ' Xa) J (x ' x ' d P (l/x; z( X )) (15) 



and 

where the superscripts denote the order of expansion in the potential. 



C^\l;za,z P ) = cf\l;za,Zp) = , (16) 



2.3. Fourier space: second order 

To work to second order, we need the usual convolution theorem for the product of two fields U and V is 

[UV](\) = [0 * y](i) = f Jj-L#(i')Pa - 1') ■ (17) 

Introducing 

9Ji(l',l-l';z s ) = fW^fV^'^ll'^^-l';/), (18) 
Jo Jo 

and using the second term from Eq. (O and the relation between convergence, rotation, shear and 0yy, the second order corrections 
to convergence, rotation and shear can be written as 

* (2) (l,z s ) = -2 r-^L[l'.(l-l')]l-l'TO(l',l-l';z s ) , o/ 2) (l,z s ) = -2 r-^L[l'.(l-l')]«'sin(^)TO(lM-l';Zs) (19) 



and 

j2i/ 



f\l, z s ) = -2 fi-L [I' • - >')] H'G/d, 10SW (1', 1 - 1'; z s ) . 

J (2w) 2 



(20) 



Here the superscript refers to the order of expansion in O, and we define G\(l, 1') = cos(0/ + 0/<) and G2(l, I') = sin(0/ + 0/*). When we 
work beyond first order in the lensing potential, the shear becomes a non-linear function of the gravitational potential <£. Hence the 
power spectrum of the shear depends on the higher order correlation functions of cD. Therefore we need the Limber approximation 
for these higher order correlation functions. For the bispectrum, Eq. ( fT4b generalizes to 

(Hh;xi)4>(h-,X2)4>(h;xi)) = (*r) 2 g„(l. + h + h) So(xi "^ d(Y ' ~* 3 W ^A-^;z0n)); (21) 
x ' X{ \Xi Xi Xi I 

and for the trispectrum, 

(0(h;Xi)<f>(h-,X2)(t>(h;XJ)<f>(h-,X4)) = (2tt) 6 D (h +1 2 +I3 +U) 2 Tq,\—, — , — , — ;z(Xi) . 

x /c X\ \Xi Xi X\ X\ I 

(22) 

where the subscript "c" denotes a connected function. 

As an example, we consider the correlation of two 9Ji functions, 

(TO (I',l - 1'; za) 9K (l'", 1" - 1"'; Zp)) = p d X j*' A X " J* d X ' J dx"'W(x,XaW(x',xW(x",XpW(x"',x") 

» (IV) ~4> (1 " IV ) 0" V) 0" - 1" V")) ■ (23) 
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The expectation value here can be broken up into a Gaussian (Wick's theorem) piece and a connected (non-Gaussian) piece. The 
connected piece vanishes because the d>D -functions in Eq. (l22l forces = x' — x" - x'" where the window functions vanish. Of the 
3 possible contractions for the Gaussian term, the only one that survives is \" — X > X'" - X' ■ Thus, 

(9K(l',l-l'u a )2R(r',l" -\"';zp)) = (2n) 4 6 D (\' + l"')6 D (l - 1' + 1" - \'")M{V, \l - \'[,z a ,Zp), (24) 

where we have introduced the mode-coupling integral 

M[l,l;z a ,Zj3)= I dx J 1 d * a — P ® ( l lx,z(x))P<s> (I IX , zix )) ■ (25) 

«7o X X 

Note that Eq. d24"l > is true even for a non-Gaussian density field. 

The third order terms each require specialized treatment, so we handle them on a case-by-case basis below. 



3. The corrections to the power spectrum 

We can now calculate the higher order contributions to the reduced shear power spectrum by Taylor expanding the reduced shear in 
terms of the shear and convergence to contain all terms up to (9(0 4 ), 

(?e/b(1,z») §e/b(I', zp)) ~ (if * (1 + k + * * K)) E/B (l,z a ) (f* (1 + % + ic* k))eib$ ,zp)) , (26) 

where * denotes a convolution, and where the shear and convergence need to be expanded in terms of the potential according to Eq. 
([8]) and projected into E/B components using Eq. 0131 . 

As the power spectra depend only on the magnitude of 1, we can choose l||x, which implies 77 ft = (1,0) and thus 7e(1) = yi (1), 
and simplifies the calculations without loss of generality. Consider for example the correction to the E-mode power spectrum arising 
from the correlation between second order corrections 

(§< 2) (i, z a ) ~gf(V,z,)) = ((r (2) + y (1) * * (1) ) E ft z a ) (r (2) + y (1) * * (1) ) E fl', zp)) 

= (rf ) (U ff )f E 2, (i',^)) + ((r ( 1 1, ^ (1) )(U») (r (1) ** (1) ) E (i'^)) 
+ (yfft za) (y (1) * * (1) ) E fl', zp)) + ({f\ l) * d, z a ) yf(\',z p )) 

= (ri 2) (D m')ff(i')) + ((f « * * (1) )a) Tvao {yf * * (1) )ao) „ + 2 (f fa) m>) (ff * *«)(i')L - ( 27 ) 



where in the last step we have rewritten the E-mode component using Eq. ( TOT l and where we define the symmetrized expectation 
value 

(A(l)B(l')) afj = i [(Aft Z a )B(Y, Zp)) + (Aft ZpW, Za))] , (28) 



to shorten our notation 



Noting yj^ft = and <y E ftyB(l')> = 0, we can expand Eq. (|26]> to <9(<D 4 ): 



(|e(1, Za) SeW, Zp)) = Cg l) (l; za, Zp) + AC^ 2) (/; za, Zp) + AC£ 3) (/; z a , zp) + AC| 2) (Z; z a , Zp) + Ac£ 3) (/; z., Z/J ) + ACg 2) ft z ff , zp 



0(0 3 ) reduced shear <3(ri> 4 ) shear 6>(U> 4 ) reduced shear 

= z ff , z,) + 2 (r»a) tvcio (rf* * a) )ci'))^ 2 (y<> wyf ao)^ (rf } ® r^ao)^ 

+ 2 {(??>(!) 77ft) (yf * ^)ft)> a , + (y^ft ^ft) (yf * * (1) )d')) ff/3 + (yfft r/d') (y? * * * (1) )a')) J 

+ {2 (yf (1) W) (yf * * (1) )(I')) a / ((??* * (1) )ft rzft) (yf* J (29) 

(isftza) § B ft,Z/0) = AC^\l;z a ,Zp) + Acf*\l; z a , Zp) 

- (rf (i) ^^a')ff d')) a/ 2 (yf ft ^0 (f^* ((r^- (1) )(i) e„r /( i') (y«* ^o)^ . 

(30) 



where we have omitted terms such as AC1 £ which vanish under the Limber approximation. 



3.1. Multiple-deflection shear corrections 

The shear-only corrections come in two flavors: the "22" (2nd order-2nd order) terms and the "13" terms. The "12" terms are 
mathematically of order <£ 3 , and hence one might expect them to be present if the matter bispectrum is non-zero. However, they 
vanish in the Limber approximation due to the W(x',x) factor in Eq. ( fT8l ). which is zero whenever^' = x- 
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The "22" B-mode shear correction can be written as 

1-1' 



{yf{\,z a )yf{\",z p )) = 2e„77(l) J ^ [1' ■ (1 - 1')] U'Gj{\,\') 2e HK T H (l") 

xJ— 2 [l'" • (i" - 1'")] i"i"'G K (\", x") (m (i,i - i'-,za) arc (r, i" - 1'"; 

= (2n) 2 6(\ + 1")4Z 2 J '^L (/' S in ^,) 2 [1' • (1 - l')] 2 M (/', |1 - \'\;z a , Zp) , (31) 



where we have used Eqs. (TOl l20l l24l i and <f>i - repeatedly. By comparison with Eq. ( fT9b one can see that AC< 22) = ACf ) 
Similarly, 



AC| 2) (/; ^) = 4l 2 J ^ (/' cos #) 2 [1' • (1 - 1')] 2 M(l', |1 - 1'|; z„, ^) , (32) 

and ^ 

ACf'(Z; = 4 J-|L(l.l') 2 [l'-(l-l')] 2 M(r,|l-l'|;z ff ,z^ . (33) 

The integrals in Eqs. ( 1321 [33} are dominated by angular scales corresponding to the peak of the matter power spectrum, which is 
at scales much larger than those typically probed by lensing: If we define l c = 1 - I', then for small / c (compared to I of lensing 
experiments) the contribution to these integrals scales as J d 2 l c / 2 cos 2 (l, l c )M(l, l c ;z a ,Zp). Assuming an effective power-law index 

n'f for the non-linear matter power spectrum P$ y „i(k), the / c -dependence of M(l,l c ;z a ,Zp) scales as l" s 4 . So the contribution to 

the integral per logarithmic range in l c scales as Z c s , which is dominated by scales corresponding to the peak of the matter power 
spectrum. 

The "13" correction in principle has three parts: those arising from the 3A, 3B, and 3C terms of Eq. ([§). Let us consider the 3B 
term first. The expectation value of the product of two Fourier modes is 



rXa rxp rxi rxi r h^i ' f AflJ' 

(^(l,ZaW?j B \L,Z P )) = 16 J d* J d Xl J dx[J d X '{ J — J ^fftoJfffc^Wt^W,*) 

xl a l h L' c L'jL' ( ;(L -V- L")«(L - V - L") C (L - V - V') d 

x(ftkxtf(L - L - L";^ 1 )0(L';y i )0(L";y i ')) . (34) 

In the Limber approximation, the only non-vanishing contraction is at^ - X\ ® a &X > \ = x'\- The (Jo-functions then enforce L' t = -L" 
and L = -1. We thus find: 



rnwnxaW rx 
b/ffi&ZaW^O+Ze)) = (2;r) 2 5 D (l + L)16 d X d X \ 

Jo Jo 



o X 2 x' 2 



x f ^ 2 P^(llx\z(j())P^{L'lx\\z(x\))lali^Wck- (35) 

The integrand is odd under L' — > -L', and hence the "13B" correction to the shear power spectrum vanishes. 

The "13C" correction is zero because the restriction^-" < x ' < X m Eq. © implies that there are no allowed contractions within 
the independent lens plane approximation. This leaves us with the "13A" correction, which is similar to "13B", except with the 
replacement L'. — > lj. The choice l||x implies that the only non-vanishing component of "13A" is (ifr^(l,z a )'pfi\l',Zp)). Hence we 
find f 

AC£ 3) (Z; za, zp) = ACg 3A) (l; z a ,z P ) = ACf A \l; z a ,z p ) = -4l 4 J (1 ■ I') 2 M(l, I'; z a , z p ) . (36) 

There is no "13" B-mode shear or rotation power spectrum because y B (1, z a ) and o< (I) (l, z a ) vanish. 

The dimensionless shear power spectrum, A?™(0 = 1(1+ l)C^\l) / (2n) 2 scales as A 2< e H) (Z) oc l n f +2 , while the corrections A?" 3 ' (/) 

and A 2 ^ 22) (0 scale as l"* B+4 . The main contribution to these corrections at large I is the bulk deflection on small scales by large 
wavelength density perturbations which causes only small local distortions. Thus the "22" a nd "13" terms largely cancel, similar 
to the perturbative calculation of the one-loop correction to the density power spectrum (e.g.|Vishniac 1983). As these corrections 
diverge for large I and have opposite sign, their numerical difference needs to be evaluated carefulljQ. 

The dotted lines in Fig. Q] illustrate their magnitude for z a = Zg — 1 using th e fittin g formula of ISmith et alj d2003l) for the 
non-linear matter power spectrum with the transfer function from Efst athiou et al.l d 19921) for the numerical integration. Here the 
combined E-mode correction is negative at small / and positive for I > 4200. These corrections are at least 4 orders of magnitude 
smaller that the linear theory result ci . 



6 Apply a variable transform 1" = 1 - 1' to AC? 2 * and cancel diverging contributions at 1" by rewriting the integral as AC? 2) + AC~' E 3) = 

4 /0 A • <J" + '» 2 0" • G" + { M ® + nr'l^Zp) - M(l,l";z a ,zp)) + ((1 • (I" + l)) 2 (1" • (1" + l)) 2 - / 4 (1 • \"f)M(l,l";z a ,Zp\ where the 
azimuthal integration of the second term can be done analytically. 
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Fig. 1. Linear order shear power spectrum (thick solid line; Eq. (fT2l i) and corrections up to <9(<E> 4 ). 

Left. The dashed/short dashed lines show the fourth order corrections to the E/B-mode shear power spe ctra that arise from re laxing 
the Born approximation and including lens-lens coupling in the calculation of the shear (Sect. [3~T1 cf. ICoorav & Hul J2002)). The 
E-mode correction is negative at small / and positive for I > 4200. The dashed-dotted line illustrates term C^ 2) (cf. Eq. ( 1321 ) which 
contributes to the E-mode shear correction , the divergency is cancelled by Eq. Q61 >. 

Right. The dashed/short dashed lines show the combined forth order corrections to the reduced shear E/B-mode power spectra 

(Sect. [3~2l Tabled- The dashed-dotted line shows the third order correcti on to the reduced shear E-mode power spectrum. 

We ass ume a source redshift z a — Zp — 1 and use the transfer function from Efstathi ou et alJd 19921). the fittin g formula of Smith et alj 
(l2003h for the non-linear matter power spectrum, and the fitting formula of Scoccimarro & Couchman ( 200 1 ) for the non-linear mat- 
ter bispectrum. This figure assumes a flat ACDM cosmology with (Q m , Qb, cr 8 , h, n) = (0.3, 0.05, 0.9, 0.7, 1) to enable comparison 
with previous calculations. 



Note that unlike the results of Co orav & Hul d2002l) . our calculations agree with the expected equivalence between the tangential 
shear and convergence (cf. Eqs. ( l32l [33l [36l l). as well as between cross shear and rotation power spectra (cf. discussion after Eqs. 
(EQEUl). 



3.2. Reduced shear corrections 

The same methodology used for the corrections to the shear power spectra can also be used to compute the reduced shear terms in 
Eq. ( 1291 . Corrections to the reduced shear power spectra which combine second order and first order distortions contribute through 
two Wick contractions, for example 

(ff(l,z a )(^*i<V)jl'',*fi)) = -2euTj{\) J^-[l'.(l-l')]U'Gj(l,l')e HK T H (l") J -^l"' 2 T K (l" -l'"\ 2 

x JVW',**) j%"'W(x"',Xfi) (hb (l', l - v-.zc) (!'";*") 4> 0" - 1"';*"')) 

= - (2tt) 2 5 d (1 + 1") 2 f ^- [I' • (1 - 1')] |I - If ll' 3 sin(^) {sin(2^) + sm(2^_y)} 

x <k 2 d x 72 p<1> " lx> z0a)P<i>(\i - 1 \lx ;zCt )) > (37) 

Jo X Jo X 

where we have used <pi — and ejjTj(l')Tj(l") = sin(20/» - 2<pi>). 

Corrections to the reduced shear power spectra which combine only first order distortions contribute through all Wick contrac- 
tions plus a connected contribution, for example 



C aIv Cttf-]'" 

((f w ^ (1) ) B (U«)(f«* * V'*>> = euTKi) J ALr /( i')^ ff (i'') fyi-Tjn 

x (^\V,z a )K m (l - l',z ff )* (1) (l"',Z/0* (I) (l" - V",zfi)) (38) 
= (2tt) 2 <5 d (1 + 1") J ^5 sin(2^){(sin(20,O + ^{l^C^l'-^Zp^iV - \;z a ,Zp) 

+ f^ S m(2^,)T K (l',l-l'X',-l-l'''-,z a ,Za,z f3 ,Zf i )}, (39) 
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Table 1. (9(<t> 3 ) and 0($> A ) contributions to the reduced shear E-mode power spectrum. 



Type 


Term 


Contribution to C g - E (/) 


0((t> 3 ) reduced shear 


$ ) Q.Za){? n ** 1) ) E <y',Zfi)) 


2 f cosf2<fc )BJl. Y.-l- Y: z„. z«. z«) cf. (Shapiro 2009) 


shear with 


(yf{Yz a )yf{Y', Zl3 )) 


4/2 / (/ ' cos ^'')) 2 d' ' d - !')) 2 M (/', |1 - 


multiple deflections 






reduced shear with 
multiple deflections 


{^\lz a )(y w *i<%a",z,)) 


-2/0(l-l') 2 / 2 /' 2 cos(20,O 




{^\lz a ){y^*i< w ) E (\",z,)) 


_ 2 f01., W | 1 + 11 




(7^fl,^)(y (1) *^ E fl^)) 


-2/0 [1' • (1 - 1')] |1 - If //' 3 cos(6r) (cos(2^) + cos(20,_ r )} 

x j^y^yr^, ™*>™**>p.(r/x; zteW*® - \'\lx'-M 


reduced shear with 






C { ^(l;z a ,Zp)al E {zp) 


single deflection 






+ J0cos(2#) fg^T&Y.Y", -l-\-V";za.zp,zi>,zid 




((f (1) *if (1, ) E (i^)(r (1) *^ 1) ) E (i"^)) 


|0 cos(20,o{(cos(2^) + cos(20 F _ l ))C< 1 E 1) (r; z^X^'dl' - l| ;Za ,^) 








+ J|j£ cos(2^or,(i', i - 1', l'", -l - r z ff , Z/J , 9 )} 


Table 2. (9(<D 4 ) contributions to the reduced shear B-mode power spectrum. 


Type 


Term 


Contribution to C g - B (/) 


shear with 






4/2 X0 ( Z ' sin (^/')) 2 V ■ d - l')) 2 M(l', |1 - l'|;^) 


multiple deflections 








reduced shear with 
multiple deflections 


(y ( i\lz a ){y w *^\(V',z P )} 


-2/0 W ■ (1 - 1')] |1 - l'l 2 //' 3 sin(^) {sin(20,O + sin(20,_,O} 

x ^ dx ^J^f d y m^A PoiVlx , z(x))PoQl _ i Vx >- z(xl 


reduced shear with 


((f (1) *^ 1) ) B az a )(r (1) *? cl) ) B (i"^)) 


/0 sin(2#){(sin(2#) + sin(2fr^)C£V;i.*)C& 1> (|r - \\;z a ,z fi ) 


single deflection 






+ J&p sin(2^)r,(l',l - V,V",-l-Y' , ;za,z a ,^,z P )} 



where we have omitted a term which only contributes to th e I — mode, and where T K (l\,l2,h, -Im', z a , z a , Zp, Zp) is the lensing 
tomography convergence trispectru m dCoorav & Hul 120011) which we model with the halo model of large scale structure (e.g., 
SeljaklHoOO; ICoorav & Shethll2002l) as summarized in Appendix [A] Here, the Gaussian contribution, which is the dominant term 
on relevant angular scales, is simply a convolution of the standard (9(0 2 ) lensing tomography cross spectra with some geometrical 
projection factors. Note that in the halo model framework the connected contribution to the B-mode spectrum is downweighted by 
the geometric projection factors, especially one-halo and (13) two-halo are strongly suppressed. The connected E-mode terms given 
in Table[T]has opposite angular symmetry and the connected part starts to dominate the signal above / ~ 8000. 

The analytic expressions for all contributions to the fourth order tangential reduced shear cross spectra are summarized in 
TableQ] Fig. Q] illustrates the numerical values of the different corrections. The fourth order reduced shear corrections of the lensing 
E-mode power spectrum reach the percent level at small angular scales and hence may be relevant for future weak lensing exper- 
iments. Reduced shear generates a small amount of B-mode power, which is about 4 magnitudes smaller than the E-mode signal, 
and is less than the level of B-mode power generated by observational systematics. 
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3.3. Relation between ellipticities and reduced shear 

The linear relation between some measure of image ellipticity and reduced shear ([TJ is only valid in the limit of very weak lensing 
(k <c 1, \y\ <K 1. In general the relation between image ellipticity and reduced shear depends on the ellipticity measure under 
consideration. As an example we consider two definitions of the complex image ellipticity here: 

s = — -e 21 * (40) 
1 + r 



and 

1-r 2 



l+r 2 



(41) 



where r < 1 is the minor to majo r axis ratio of the image , and <p is the position angle of the major axis. The latter is frequently em- 
ployed in observational studies dBernstein & Jarv is 2002), the former is more of theoretical interest due to its simple transformation 
properties. The full relation between ellipticity and complex reduced shear g = gi + igi is given by 

e (s) + g + 2g + g 2 e (s - ) * 

and e= - , „ * * ,^ s , (42) 



1+g'eGO 1 +|g| 2 +2^(ge< s )')' 

where 7?(z) is the real part of a complex number z, e (s) and e (s) are the intrinsic ellipticiti es of the source and wher e we only consider 
\y\ < 1, which is certainly true for cosmic shear. The linear relation (e) = g is exact (Seitz & Schn eider]| 19971). as can be sh own 
using the residue theorem. In the second case, using a Taylor expansion (Schne ider & Seitz! 19951; iMandelba um et al.ll2006l) . the 
ellipticities can be written as 

<e> = c,g + c 3 |g| 2 g + 0(g 5 ) * (2 - e (s)2 ) g + (-2 + 5e (s)2 - 3e (s)4 ) |g| 2 g, (43) 

where e (s) is the absolute value of the intrinsic ellipticity of the source galaxies. In the practical case of a distribution of intrinsic 
source ellipticities, one should replace the powers of e (s) by their moments (e (s) "). Shear is typically estimated by taking the mean 
observed ellipticity (e) and dividing by the response factor c\. To (9(0 4 ), this shear estimator reads 

g= — =g+-|g| 2 g. (44) 

Cl Cl 

The last term gives rise to one additional contribution to the power spectrum of g e: 

2^ Off&Za) (f (1) * 7 (1) * f (1) L(l',^)) =1- SuTS') f^- 2 f ^ 2 T„(l'")T H (l")TAV - 1" - 1"') 

x (^«(i, Ztt )^a'",z^ (1) a",z^ (1) a' - \" - r,^)) 

= (2tt) 2 6 d (1 + 1') 2^ J ^l{(2cos 2 (20 r ) + \)C^\l,z a ,Zp)C^\r-z P ,Zp) 

+ J ^ cos(2^ - 2<p lt „) cos(2^_ 1 _ F _ lw )n(l, 1", 1"', -1 - 1" - 1"'; z a ,Zp, Zp, z P )) 

= (2jt) 2 6 d (1 + V) 2 C M [ 2C^\l;z a ,zp)cr\ E (z P ) + J^^- cos(20,„ - 2^)co S (20_i_ 1 »_ 1 »)n(l,l",l'",-l - 1" - V" ;z a ,zp,zp,zp)), 

(45) 

where we have performed the angular integration of the Gaussian contribution in the last step and introduced the shear dispersion 

4M = f^rc£\i'-,ze,z P ). (46) 



For the case of the s ellipticity, linearity implies c\ — 1 and C3 = 0. In this case, the correction of Eq. d45l > vanishes. For the case of 
the e ellipticity, we have 

c 3 -2 + 5e (s)2 - 3e (s)4 

7, - 2 -et* ■ (4?) 
The magnitude of this corrections for the e ellipticity with (e (s)2 )^ 2 = 0.6 is illustrated in Fig. [2] 
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Fig. 2. Linear order shear power spectrum (thick solid line; Eq. j!2b ) 
and (9(0 4 ) lensing bias and ellipticity estimator corrections. 
The short dashed (dashed) lines show the lensing bias corrections to 
the B-mode shear power spectrum (Eq,l60t assuming q=l (q=2). The 
dotted (dashed-dotted) lines show the lensing bias corrections to the 
E-mode shear power spectrum (Eq.|59l> assuming q=l (q=2). 
The fine solid line illustrates the magnitude of the correction arising 
from the conversion between ellipticity and reduced shear Eq. J45b for 
the e ellipticity with (e <s)2 ) 1/2 = 0.6. This correction is negative and 
its normalization depends on the distribution of source galaxies (see 
Sect.[33]for details). 

This figure uses the same cosmology and source redshifts as Fig.Q] 



3.4. Lensing bias corrections 

Galaxies are only selected for shear measurement if they are large enough and bright enough to measure their shape. As lens- 
ing changes the observed brightness and size of the lensed galaxies, the number of galaxies selected above some magnitude 
and size threshold is correlated with the lensing field (this is the well know magnification bias, and the size bias discussed in 
dSchmidt et al.ll2009bl)). Neglecting source c lustering, the normalized observed galaxy overdensity due to lensing magnification is 



(scnmiatet al. ZUUVb)). iNeglecting source cms 
given bv dSchmidt et alJl200alHui et alJl2007h 



1 + qkCo) + C\k z (n) + Cz\y\ (n) / , i n \\ i , i ~>w 

1 + <Wn) = \ +Cl {#)Vc 2 (\y?) * 1 + qK(n) + Cl ^ (n) ~ + Cl (D) ~ < W » ' (48) 

where we expanded the magnification to second ordefl, and where C\ = q(q + l)/2 and C 2 = q/2. The parameter q is de termined 
by the slope of the luminosity and radius distribution of the sample galaxies and typically q ~ 1 — 2 dSchmidt et al.l l2009b). 

Hence the sampling of the shear field measured from galaxy pairs is modulated by the lensing magnification implying that the 
observed shear depends on the true shear and the galaxy overdensity 

g 7 0bs (n) = g,(n) (1 + «Wn)) * yi(n) (l + <n) + * 2 (n) + q K (n) + ^ 2 (n) + d[K 2 (n) - <* 2 >] + C 2 [\y(n)\ 2 - <| r | 2 )]} . (49) 

The standard pair b ased estimator for the reduced shear correlation functions £; a t — {g a gb) then becomes (for details see 
ISchmidt et ai1l2009ah 

{hm) = I ^gi bs (n)gf s (n + 8)\ , (50) 

where N is the observed number o f galaxy pairs with separ ation 9 relative to that expected for a random distribution; this is just the 
^ correlation function estimator (iPeebles & Hauserll 19741) . For large-angle surveys, TV converges to the correlation function, 

N -» 1 + <<5ien S (n)(Wn + S)). (51) 

Therefore we may write 

,„ . (gf> s (n)gf s (n + 6)) 

V ' 1 + <5len S (n)5lens(n + 6)) 

This can be converted to products of correlation functions by conversion to a geometric series, 

oo 

(hM) = {g? s (n)gf s (n + 6)) 2(~D y <4* s (n)<Wn + 8))"; (53) 

v=0 

we then note that the v term in this expansion is of order 0{<b 2+lv ). Since (|/y(0)) is desired to <9(<t> 4 ), it suffices to keep only 
the v — and v — 1 terms. Moreover, in the v — 1 term, we only require the lowest-order expansion of the correlation function 
<<5iens(n)<Si ens (n + 0)), i.e. 

<<Wn)(5i ens (n + 0)> = q 2 ( K (nMn + 0)} + 0(<D 3 ). (54) 

We also need only the lowest-order expansion of ^g° bs (n)g° bs (n + 6*)^ in the v — 1 term, i.e. we can approximate it as {jj(n)yj(n+8)). 
Thus we reduce Eq. d52b to 

(hM) * (sf s (n)sf s (n + 0)) - q 2 <r/(n)ry(n + 0)) «nMn + 8)) . (55) 



7 We note that (k) = 0. This is because by rotational symmetry the mean deflection angle (d) = 0, and therefore its derivative Wi t ) = 0. 
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A straightforward generalization to cross-correlations between different redshift slices gives 

(iu(0, za, Zfi)) * (s° bs (n, z a ) 8 f s (n + 6, zp)) - q 2 (r/(n, z a )yj(n + 0, zp)) (k(u, z a )i<(n + 6, zp)) . (56) 

We now turn to practical computation. The terms involving ^g° bs (n,z (y )g° bs (n + 0,zp)) are identical to terms that we have 

calculated previously, except with additional factors of q, q 2 , C\, and/or C2, and hence present no new difficulties. The final sub- 
traction term is the product of two expectation values and hence is different from terms that we have previously considered. This 
"product correction" can be evaluated by noting that its contribution to the observed correlation function is the product of the shear 
and convergence correlation functions. In Fourier space, this means that its contribution to the power spectrum is the convolution of 
the shear and convergence power spectra: 



AC^(l) = -q 2 J -0yC ym (l')C K Ml - I'D, (57) 



where all power spectra carry the redshift indices z a , Zp- Specializing to the case where 1 is along the x coordinate axis, and recalling 
that the /s-mode shear and convergence power spectra are equal, we can then infer a contribution to the observed E'-mode power 
spectrum 

AC£ od (l) ^-q 2 j^L cos 2 (2« C™(Ocg x \\l - 1'|); (58) 

the B-mode contribution is similar except for the replacement cos 2 — > sin 2 . 

Similar to Eqs. (l26l l29l . we now expand (g E /Bd)g E / B d'')) to ^ nc ^ t ^ le f° urm order power spectrum corrections ACtj 8 ^ which 
arise from lensing bias 



ACif (/; z a ,zp) = (2^(1 + 1") 2q f(^>(l, z a ) (y m * * (2) ) E d", zp)) + (y^d, z a ) (r (2) * * (1) ) E d", zp)) + (y ( i\l z a ) (r (I) * * (1) ) E d"' zp))] 

+ (2q + q 2 ) ~^) E (i Za ) (r (I) * * (1) ) E d"^) + 2C 2 [(y£\iz a ) (r (1) * y (1) * r (1) ) E d",^)) c + c ( } E l \i-z a ,zp)cr 2 E (zp)] 
+ 2q(yV\i Zo ) (y (I) * * (1) ) E d",z,)) + 2C, (y«\i,z a ) (y (I) * ^ 1) ^ (1) ) E d",^))J 

C ' a/3 

and 

AC^(l;z a ,Zp) = (2n) 2 6 D (l + I")( 2q(yf{\,z a ) (y (1) * k {1 \{\" ,zp)) + (2q + q 2 ) ((y (l) * K ( \{l,z a ) (y (I) * k { \(1" ,zp))\ 

\ " * " lap 

~ ql I S sin2 ( 2 ^)0 /,;z -^ )C ^ 1)(|1 ' - • (60) 



In Eq. (1591 1 we have simplified the terms which involve the variance of shear or convergence, e.g. the term in Eq. (1561 ) which is 
proportional to C\ becomes 

Ci (yV&Za) (y (I) * * (1) ** (1) ) E d'^)) - C, (y^(l,z a )y^(l',zp)) (fty) = C, (y£\l,z a ) (y (I) * * (1) ) E d») c • (61) 

Here the second term is canceled by the disconnected part of the first first term arising from the Wick contraction C\ (yy) (kk), the 
two other Wick contractions of this term vanish after azimuthal integration. An explicit expression for the connected term is given 
inTableQ] 

For the redshift range and cosmology considered in this paper, the second term and third in Eq. (f60l > are the dominant contri- 
butions. These terms partily cancel and on scales / > 50 lens ing bias effectively in creases the B-mode power spectrum by approxi- 
mately a factor ( 1 + 2q), which is smaller than the findings of Schmi dTet al.l (l2009al) who only considered the Gaussian contribution 
to the second term in Eq. d60b . T he B-mode sign al is largest for small angular scaled and high source redshifts. Assuming q < 2 
and a WMAP5 cosmology (Kom atsu et alJl2009l) . for sources at z < 3 and in the range I < 10000 the B-mode power spectrum is 
suppressed by at least a factor 500 (a factor 3000 for z < 1) compared to the shear E-mode po wer spectrum. 

Lensing bias gives rise to a third order correction discussed by Sch midt et al.1 (l2009al) . which is q times the reduced shear 
correction analyzed by Shapiro (2009). The fourth order E-mode correction generated by lensing bias Eq. d59b is more complicated 
and we will discuss its impact on the E-mode power spectrum in Sect. [4] 

The lensing bias E-mode and B-mode corrections are illustrated in Fig. [2] assuming a source redshift z a = Zp = 1. Due to 
uncertainties in modeling the non-linear clustering of matter on small scales we restrict our analysis to / < 3000, on these scales the 
lensing bias corrections are below 1%. 
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Table 3. Z values for the (9(0 4 ) corrections for different ellipticity estimators with lensing bias. 



estimator 


q = 


q=l 


q = 2 


e 


1.14 


3.19 


5.31 


e, <e^>" 2 = 0.6 


0.12 


2.13 


4.25 



4. Impact on future surveys 

The corrections derived in Sect. [3] generate a small amount of B-mode power, and have a < 1% effect on the ellipticity E-mode 
power spectrum. These are well below the error bars of current surveys and therefore have no significant effect on published results. 
However, future "Stage IV" surveys such as LSST, JDEM, and Euclid will be sensitive to sub-percent effects. We can quantify 
the importance of the higher order lensing corrections by comparing the corrections to the power spectrum AC(/;z ff ,z l g) to their 
covariance matrix. Quantitatively, 



Z= Yj iCov- 1 [C(l;z a ,z^),C(l';z a ',z fi ')]}AC(l;z a ,Zp)AC(l';z a ',z^) (62) 

represents the number of sigmas at which the corrected and uncorrected power spectra could be distinguished by that survey. 
Corrections with Z <k 1 are negligible in comparison with statistical errors, whereas corrections with Z » 1 must be known to 
high accuracy to make full use of the data set. We have computed Eq. d62l assuming a WMAP5 cosmology dKomatsu et al.1 12009) 
for a model survey with a surface density of 30 galaxies/arcmin 2 , median redshift z me d = 1.1, and sky coverage of 10 4 deg 2 , as 
appropriate for some of the proposed versions of JDEM. The power spectra were computed in 14 redshift slices and 12 /-bins with 
a maximum multipole of / max = 3000. The algorithm fo r computing the cova riance matrix is as described in Appendix A.2.d of the 
JDEM Figure of Merit Science Working Group report (lAlbrecht et al.ll20 09). Without lensing bias (q = 0), we find Z = 1.14 for 
the linear ellipticity estimator e; for the standard estimator e and for an rms ellipticity0<e <s)2 > 1/2 = 0.6, we find Z = 0.12. Including 
the lensing bias corrections from Sect. 13. 41 increases the significance of the corrections as detailed in Table [3] Note that the table 
includes only the Of® 4 ) c orrections, and does not include the (9(<t> 3 ) corrections that have previously been considered ( Shapiro 
2009; Schmidt et al. 2008). Thus, the perturbative corrections to the weak lensing approximation are expected to be at the level of 
~ 1 - Act. These corrections will have to be taken into account for future surveys, but given that they are only ~ 1 - 5cr and should be 
accurately calculable (either directly via ray-tracing simulations, or by analytic expression in terms of the moments of the density 
field, which can be determined from A-body simulations), they should not represent a fundamental difficulty. 



5. Discussion 

We have calculated the reduced shear power spectra perturbatively to fourth order in the gravitational potential, accounting for the 
differences between shear and reduced shear, relaxing the Born approximation, and including lens-lens coupling in the calculation 
of shear and convergence. The full set of corrections to the reduced shear power spectra are given in Table [TJ (E-mode) and Table 
[2] (B-mode). The ellipticity power spectrum contains additional contributions, Eq. (1451 . which arises from the non-linearity of the 
shear estimator and depends on the specific definition of ellipticity used, and Eq. (l59l which is caused by lensing bias. Through 
order O 4 , this is the full set of corrections to the power spectrum arising from the lensing process itself. All corrections have been 
derived within the Limber approximation, and the analysis of "12" type multiple-deflection corrections is left for future work. Other 
corrections associated with the source galaxy population, such as source clustering and intrinsic alignments, are not treated in this 
paper. We find that, depending on the properties of the source galaxy population and on the type of shear estimator used, these 
corrections will be at the ~ 1 - 5cr level, and thus should be included in the analysis of future precision cosmology weak lensing 
experiments. 

That said, we caution that there are other areas in which the theory of weak lensing needs work if it is to meet ambitious 
future goals. Current fitt ing formula of the non-linear dark matter power s pectrum have an ac curacy of about 10% at arcminute 
scales (ISmith et al"1l2003l) and the uncertainty exceeds 30% for / > 10000 (iHiibert et aT1l2009l) . due to this difficulty in modeling 
the non-linear gravitational clustering angular scales of I > 3000 are likely to be excluded from parameter fits to cosmic shear 
measurements. Utilizing near-future A-bo dy simulations it will becom e possible to determine the non-linear dark matter power 
spectrum with percent level accuracy (e.g.. lHeitmann et aDl2008l 120 09). However, this does not account for the effect of baryons, 
which will likely be important at halo scales and depend critically on the details of baryonic processes (cooling, feedback) involved. 
Baryons in dark matter halos which are ab le to cool modify the structure of the dark matter halo through adiabatic contraction 
(Blument hal et al.lll986UGnedin et afl l2004). ca using deviations of the inner halo pro file from the simple NFW form and changing 
the halo mass - halo concentration relatio n (e.g. lRudd et al . 2008; P edrosa et al.ll2009l) . The latter can be constrained though galaxy- 
galaxy lensing (Mandelbaum et al. | 2006|). or could be internally self-calibrated in a weak lensing survey via its preferential effect on 
the small-scale power spectrum ( Zentner et al.ll20 08). Baryons in the intergalactic medium may make up about 10% of the mass in 
the universe, and if their distribution on Mpc scales has been stron gly affected by non-grav itational processes then they could pose 
a problem for precise calculation of the matter power spectrum (see lLevine & Gnedinl20 06. for an extreme and probably unrealistic 
example). 



The rms ellipticity here includes both the + and x components, so it is V2 times the rms per axis. 
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Given these uncertainties in modeling the non-linear matter distribution and that all the corrections derived in this paper are 
integrals over the non-linear matter power spectrum, bispectrum and trispectrum, we refrain from calculating (9(<D 5 ) and higher 
corrections. We expect that the corrections derived in this paper are sufficient to model the perturbative relation between the non- 
linear matter distribution and the lensing distortion in weak lensing surveys for the forseeable future. 
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Appendix A: Halo Model Trispectrum 

The trispectrum T(k] , k2, k3, kzt) of the dark matter density contrast 6(k) is defined as 




(A.l) 
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We model the dark matter trispectrum using the halo approach (Seljak 2000t ICoorav & Sheth 2002), which assumes that all matter 
is bound in virialized structures, which are assumed to be biased tracers of the density field. Then the statistics of the density field 
can be described by the dark matter distribution within halos on small scales, and is dominated by the clustering properties of halos 
and their abundance on large scales. In this model, the trispectrum splits into four terms, which describe the 4-point correlation 
within one halo (the one-halo term r lh ), and between 2 to 4 halos {two-, three-, four-halo term) 

t = T lh + (r 2 2 2 h + r 2 3 h ) + r 3h + r 4h . (A.2) 

The two-halo term is split into two parts, representing correlations between two or three points in the first halo and two or one point 
in the second halo. 

As halos are the building blocks of the density field in the halo approach, we need to choose models for their internal structure, 
abundance and clustering in order to build a model for the trispectrum . In the following w e summarize the main ingredients of our 
implementation of the halo model convergence trispectrum following (jCo orav & Hu 200JJ). 

We assume the halo profiles to follow the NFW profile dNavarro et al.l l 1997) 

, s = A vir pc 2 1 

Ar ' C> 3(ln(l + c) - c/(l + c)) r/r m (1 + cr/r vir ) 2 ' 

where A v [ r and p are the d ensity contras t and m ean density of the universe at virilization, and c(M, z) i s the halo concen t ration , 
which we model using the Bullo ck et alj d2001l) fitting formula. We model the halo abundance using the She th & Tormenl fl999) 
mass function 

— dM = £-f(v)dv = —A fl + (avy p ] V^expf-— ) — , (A.4) 
dM M M v \ 2 1 v v 

where A and p are fit parameters, and v is the peak height v = S c /(D(z)cr(M)). cr(M) is the rms fluctuation of the present day matter 
density smoothed over a scale R = (3M/4^p)'^ 3 , and D(z) is the growth factor. To describe the biased relation between the dark 
matter halo distribution and the density field, we assume a scale independent bias and use the fitting formula of ISheth & Tormenl 
(1999) 

Kv) = 1 + ^ + - , (A.5) 

5 C 5 c (\+{av)P) 

and neglect higher order bias functions (b 2 , etc.). Following the notation of ICoorav & Hul d2001l) we introduce 

,^;z) = f dM ^(j) b l3 {M)p{h,c{M,z))---p{k tl ,c{M,z)), (A.6) 

which describes the correlation of p points within the same halo, and where bo = 1 and b\ is given by dA.5b . Then 

T lh (k l ,k 2 ,k i ,k 4 ;z) = I° 4 (k l ,k 2 ,k 3 ,k 4 ;z) (A.7) 

T^(k u k 2 , k 3 , k 4 ; z) = PUkx)D{z)ll(k 2 , k 3 ,k 4 ; z)l\(k x ;z) + 3 perm. (A.8) 

r 2 2 2 h (ki,k 2 ,k 3 ,k4;z) = P^(k n )D(z)l\{h,k 2 -z)l\(h,k 4 -z) + 2 perm. (A.9) 

r 4h (k 1 ,k 2 ,k 3 ,k 4 ;z) = r pt (ki,k 2 ,k3,k 4 ;z)/ 1 1 (fei;z)---/ 1 1 (fe4;z) , (A. 10) 

where k fl £ = k fl + k/,. We neglect the 3-halo term, as it has neg ligible effect on our calculation, and simplify the 4-halo term using 
just the trispectrum given by perturbation theory T pt dFry||1984h . 

Finally the tomographic convergence trispectrum can be written as 

2222 f W(x,Xa)W(x>XBW(x>XrW(x,Xs) 
z 7 , z s ) = l\l 2 ltl m d X — b : FoGi IX, h/X, hlX, -I123 lx\ zCf)) 

= \-a m H^\ Jd x — 6 {l+z(x)) T(\i/x,h/Xrh/X,-h23/xiz(x)), (A.ll) 

where we have used the Poisson equation to relate the potential trispectrum to the matter density trispectrum. 



